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We analyze the sensitivity of a network of interferometer gravitational-wave de- 
tectors to the gravitational-wave stochastic background, and derive the dependence 
of this sensitivity on the orientations of the detector arms. We build on and extend 
the recent work of Christensen, but our conclusion for the optimal choice of orienta- 
tions of a pair of detectors differs from his. For a pair of detectors (such as LIGO) 
that subtends an angle at the center of the earth of ^ 70° , we find that the optimal 
configuration is for each detector to have its arms make an angle of 45° (modulo 90°) 
with the arc of the great circle that joins them. For detectors that are farther sepa- 
rated, each detector should instead have one arm aligned with this arc. We show that 
the broadband sensitivity to the stochastic background of a detector pair which are 
^ 3000 km apart is essentially determined by their relative rotation. Their average 
rotation with respect to the arc joining them is unimportant. We also describe in de- 
tail the optimal data-analysis algorithm for searching for the stochastic background 
with a detector network, which is implicit in earlier work of Michelson. The LIGO 
pair of detectors will be separated by ~ 3000 km. The minimum detectable stochastic 
energy-density for these detectors with their currently planned orientations is ~ 3% 
greater than what it would be if the orientations were optimal, and ~ 4 times what 
it would be if their separation were ^ a few kilometers. [The detectors are chosen 
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to be far apart so that their sources of noise will be uncorrelated, and in order to 
improve the angular resolution of the determinations of positions of burst sources.] 
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I. INTRODUCTION AND SUMMARY 



A. Background and motivation 



Construction will begin soon on the American Laser Interferometer Gravitational Wave 
Observatory (LIGO) lU, and on its French/Italian counterpart, VIRGO 0. Early in the 
next century there will likely be in operation a worldwide network of detectors, with sites 
in America, Europe, and possibly Japan and Australia [Q. It is important at this stage 
for physicists to look ahead and identify the types of science that the community might 
focus on using this network when it reaches a mature stage, perhaps a decade after the first 
gravitational waves are detected. One of the reasons for doing so is that some properties 
and parameters of the network, which ultimately will constrain what it can accomplish in 
the future, are being finalized today. The orientation of the detector arms is one example. 

One of the long term aims of this detector network will be to place upper limits on (or 
perhaps detect) the energy density of a stochastic background (SB) of gravitational waves. 
This background would be analogous to the relic 3°K electromagnetic background, except 
that its spectrum is not expected to be thermal. The spectrum is usually characterized 
by a quantity ^g{f) which is the gravitational- wave energy density per unit logarithmic 
frequency, divided by the critical energy density pc to close the Universe: 



Some possible sources of a SB include: (i) random superposition of many weak signals from 
binary-star systems Q], (ii) decaying cosmic strings |^ and first order phase transitions 
in the early Universe, and (iii) parametric amplification of quantum mechanical zero-point 
fluctuations in the metric tensor during inflation See Refs. [|IO],|lT],|l2| for reviews. The 

predicted wave strengths from all of these stochastic sources are highly uncertain, reflecting 
our relative ignorance of the relevant physics and/or astrophysics. Hence, detecting or 
placing upper limits on the SB can bring us valuable information, particularly about the 
very early Universe. 



1 dE 
Pc d?x (i(ln /) 




3 



Relic gravitational waves produced during inflation are particularly interesting, because, 
as Grischuk has shown ||13[, the energy spectrum for these waves contains a unique imprint 
of the time evolution of the Universe's scale factor a{t). We now discuss what is known 
about the magnitude of the contribution to ^g{f) from these waves, at frequencies relevant 
to LIGO/VIRGO. The predictions for Qg from cosmological models are not very flrm: they 
can vary between ~ 1 and ~ 10~^^ or less. However, observational upper bounds on ^g{f) 
in various frequency bands give interesting constraints on the models P]. In turn these 
constraints can be used, within the context of particular inflation models, to place upper 
bounds on the contribution of relic gravitons to the value of ^g{f) in the LlGO/VlRGO 



waveband, see, e.g., Ref. |TJ]. The strongest such upper bound comes from matching the 
normalization of the scalar and tensor fluctuations produced during inflation to the recent 
COBE measurement of the microwave background anisotropy [p|,p!5|. The result of this 
matching is somewhat discouraging: standard exponential inflation models predict that 
^gif) < 3 X 10-1^ at the 95% confidence level 0, far too small to be detected [cf. Eq. m 
below]. While it is far from certain that exponential inflation is correct, it seems unlikely 
that the expansion during the inflationary era was so much faster than exponential as to tilt 
the gravitational wave spectrum enough to give a detectable signal at high, LIGO/VIRGO 
frequencies. It is also possible, of course, that the observed microwave anisotropy was caused 
by physical processes of structure formation other than inflation. 

Despite these pessimistic prospects for the detection of relic gravitons by LIGO/VIRGO, 
it is certainly possible that there will be a detectable signal from other sources such as cosmic 
strings 0. Hence, it is important to determine how the detector arm orientations, which will 
not be changeable in the future, affect the sensitivity of the detector network to the SB. To 
do so is the flrst of the two principal purposes of this paper. This issue was flrst considered 



by Michelson |T6[ , and has been extensively discussed by Christensen |T^JT3l • Essentially we 
build on and extend slightly their analyses. Our conclusions are also slightly different from 
those of Christensen. 

The second principal purpose of this paper is to spell out the optimal data-processing 



procedure for searching for the SB with a network of detectors. The algorithm for two 
detectors is imphcit in Michelson [|I6| (and is incorrectly treated in Ref. we give a more 
detailed description and a generalization to a network of detectors, taking into account the 
possible effect of correlated sources of noise. We now turn to a description of our results 
and an explanation of how they relate to earlier work. 



B. Detection of the stochastic background 

The effect of the SB on a gravitational-wave detector is essentially to produce a small 
contribution to the random, gaussian noise in its output. For one detector this contribution 
will be swamped by the detector's own sources of noise, unless the SB strength is implausibly 
large {Qg ~ 10~^; see Sec. H). For two detectors which have no common sources of noise, 
however, the only contribution to the correlated fluctuations in the detector outputs will 
be the SB. By cross correlating the outputs of the detectors, the SB can in principle be 
measured. If one had identical, LIGO type detectors at the same site, oriented in same way 
so that they respond in exactly the same way to the SB, and with levels of intrinsic noise 
corresponding to the "advanced detectors" of Ref. [jll, then cross correlating would give a 
sensitivity to Qg of the order of 10"^° in the frequency band lOHz ^ / ^ 1000 Hz [1T0| . 

For a pair of separated, non-aligned detectors, two new physical effects complicate the 
analysis [|12],|16|,|171 . First, if the detectors are not aligned the same way, they will respond to 
different polarization components of the SB. Orthogonal polarization components of the SB 
are expected to be statistically independent, and so the cross correlation will be reduced. 
Second, for each mode there will be a time lag between exciting the first detector and the 
second detector, and hence phase lags in the cross correlation. For the LIGO detectors 
separated by ~ 3000 km and having maximum sensitivity at a frequency of / ~ 70 Hz, a 
typical phase lag will be of order unity. Hence, there will be some destructive interference 
between parts of the cross-correlation that are due to modes which propagate in different 
directions. Thus we expect a reduction in the sensitivity of the detector pair to the SB. 
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To analyze this reduction, it is necessary to (i) determine how to optimally process the 
output from the detectors, and (ii) find the signal to noise ratio (SNR) that results from 
this method of filtering. We also want to (iii) determine how the optimal SNR depends on 
the detector orientations, and (iv) find those orientations that maximize the SNR. Steps (i) 
and (ii) were analyzed in Refs. |T^,[T^. They found that, when optimal signal processing is 



used, the square of the signal to noise ratio for a broadband measurement of the SB is |]18 

Here f is the duration of the measurement, and S'„(/) is the spectral noise density in either 
detector. The key quantity appearing in this equation is the dimensionless function 7(/), 
which we call the overlap reduction function. It characterizes the reduction in sensitivity to 
the SB of the detector pair at frequency / that is due to their separation and non-optimal 
orientations, and its value is unity for coincident, aligned detectors. In Refs. []TB|,|T^ a for- 
mula for the overlap reduction function was derived, which expresses it as an integral over 
all solid angles of the complex phase lag between the detectors, weighted by combinations 
of the detector beam pattern functions [cf. Eq. (|2.13|) below]. Christensen [0,0 numeri- 
cally calculated this function for various detector configurations and discussed some of its 
properties. However its dependence on the detector orientations was not apparent. 

In this paper we derive an analytic formula for the overlap reduction function. Using this 
formula we are able to carry through steps (iii) and (iv) outlined above. We also determine 
how good are the choices that have been made for the orientations of the detectors in LIGO, 
VIRGO and GEO (an as-yet-unfunded British/German collaboration); i.e., we determine 
how their sensitivity to the SB compares to the sensitivity they would have if they were 
optimally oriented. 



C. Effect of detector orientations 



Our results are as follows. Call ai the angle between the bisector of the arms of the 
first detector and the arc of the great circle that joins the detectors, and similarly define (T2 
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for the second detector. Let S = {(Ti — o"2)/2 and A = (cii + 0-2) /2, so that S describes the 
relative rotation of the detector pair, and A describes their average rotation with respect to 
the hne joining them. Then in Sec. |V] below we show that the optimum SNR ( |1.2| ) is given 
by 

^2 

— = A 008^46) + 2B cos(4(5) cos(4A) 

+ Ccos2(4A), (1.3) 

where the quantities A, B, and C are independent of 6 and A, and A and C are positive. 
There are thus two possibilities for the optimum detector orientation, depending on the sign 
of B: (I) cos(45) = — cos(4A) = ±1, corresponding to each detector having an arm along 
the line joining them, and (II) cos(45) = cos(4A) = ±1, corresponding to each detector arm 
being at an angle of 45° (mod 90°) to this line. In Fig. |l| we plot the SNR for both of these 
choices of orientation, as a function of the angle /? subtended between the detectors at the 
center of the earth. The configuration II is optimal for (3 ^ 70°, while configuration I is 
optimal for detectors which are further apart. Fig. |I| shows that detectors which are close 
together are the most sensitive; the sensitivity of LIGO is roughly ~ 25% of what it would 
be if its detectors were coincident. It also shows that detectors whose planes are roughly 
perpendicular (/? ~ 90°) have poor sensitivity, as we would expect. 

We also show in Sec. that the narrowband sensitivity of the detector pair near a given 
frequency /, which is proportional to |7(/)|, is also always optimized at either configuration 
I or II. For example, at very low frequencies, |7(/)| becomes essentially the overlap of the 
polarization tensors of the two detectors [cf. Eq. ( pSj) below with P2(0) = P3(0) = 0], which is 
maximized in configuration I. This low frequency limit was previously derived by Christensen 
p^ . Motivated by this, he suggested that configuration I was always the best orientation 
to choose. Fig. |1| shows that, though this is not true for some values of the separation angle 
(3, the amount lost by choosing I rather than II is never more than a few percent in SNR. 

We now discuss the orientations that have been chosen for the detector systems that 
are under construction or that have been proposed. The relative rotation angle S for a pair 



of detectors essentially determines whether the detectors respond to different polarization 
components, or to the same polarization component, of the gravitational wave field. The 
advantage in responding to different components {26 ~ 45°) is that more information can 
be extracted from incoming burst signals. On the other hand, if the detectors respond to 
the same component {6 ~ 0°) then the detection signal to noise threshold is reduced, i.e. the 
fact that the same waveform is seen in both detectors means that one can be more confident 
that a candidate event is not due to detector noise. These considerations guided the choices 
of the presently planned values of 6 for the detector pairs LIGO/LIGO (there will be two 
LIGO detectors) and VIRGO/GEO, which are 6^0° and 26 ^ 45° respectively 

Within the context of these constraints, a key issue that we wanted to understand was 
the following: given the above values of 6, how much does the broadband sensitivity of 
the detector pair depend on A, i.e. by how much can the SNR be reduced if A is chosen 
arbitrarily instead of being optimally chosen. The answer we obtain [Sees. |^ and below] 
is that for relatively close detectors with (3 ^ 30°, the dependence on A is very weak; but 
the dependence is strong for detectors on different continents. Hence for LIGO's parallel 
detectors, the sensitivity will be close to optimal irrespective of the value of A; the present 
choice of A = 28.2° implies that the SNR for LIGO is ~ 97% of the sensitivity at optimal 
orientation. The VIRGO/GEO orthogonal detector pair, however, will have a SNR of less 
than 10~^ times the optimal value, irrespective of the value of A, if 26 is chosen to be 45° 
(as is planned, so as to optimize the information obtainable from burst sources). 

Finally, we estimate the 90% upper confidence limit that can be placed on ^g{f) by the 
so called "advanced detectors" in LIGO in a broadband measurement using a one third 
of a year integration time, to be ~ 5 x 10^^° in the frequency band 20 Hz ^ / ^ 70 Hz 
[cf. Eq. ( |6.5| ) below]. This is a little worse than earlier estimates which assume that the 



detectors are coincident and aligned [|T0,O]. 
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D. Organization of this paper 



The layout of this paper is as follows. In Sec. we define the cross correlation matrix 
for a network of detectors and give the formula for the contribution to this quantity from 
the SB. In Sec. ^ we describe the general optimal data-processing strategy, the justification 
of which is given in Appendix A. In this appendix we also derive the generalization of the 
signal to noise formula for optimal signal processing ( |1.2| ), to a network of more than two 
detectors, but more importantly to a network that has more than one interferometer per 
site: we show how to take into account the effect of correlated noise in detectors at the same 
site by introducing the concept of the effective spectral noise density of a detector site. In 
Appendix B we derive the formula for the overlap reduction function, and we describe some 



of its properties in Sec. IV. 



Next, in Sec. 0, we show how to optimize the orientations of a pair of detectors, both for 



narrowband and for broadband measurements of the SB. Section ^ describes the implica- 
tions of our results for the LIGO, VIRGO and GEO detector facilities. Finally in Sec. |VI1 
we summarize our main results. 

We use units throughout in which the speed of light c and Newton's gravitational constant 
G are unity. 

II. DETECTOR CROSS CORRELATION MATRIX 

The effect of a stochastic background on a detector network will be essentially to produce 
statistical correlations between the outputs of the various detectors. A key result which we 
will need is an expression for these correlations in terms of the spectrum ^g{f) of the 



gravitational background. This was first given by Christensen [12,17|, although it is implicit 



in the work of Michelson . We briefly describe the derivation in this section, and we lay 



the foundations for our analysis of Appendix A and Sec. Ill 



A detector network with detectors will have outputs 
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ha{t) = K^--\t) + na{t) + Sa{t), (2.1) 

for 1 < a < A^. Here ha is the strain amplitude that we read out from the ath detector; 
it consists of an intrinsic detector noise Ua, a contribution from the SB and possibly 
a contribution h^^^r^s^ from non-stochastic gravitational waves (bursts and periodic waves). 
The noise Ua and the SB induced strain Sa are independent random processes, which we 
assume to be gaussian and stationary. 

The detector correlations can be described by the correlation matrix 

Ch{r)ab = {ha{t + T)hit)) - {Kit + r)) (hit)), (2.2) 

where (...) means an ensemble average or a time average. The Fourier transform of the 
correlation matrix, multiplied by two, is the power spectral density matrix: 

/oo 
dTe''''f^Cn{T)ab. (2.3) 
-oo 

This is a positive definite hermitean matrix which satisfies the equations 

{KifYhur) = l^if - r)SuU)ab, (2.4) 

and 

(^^i j dtw,{t)K{t)\^ = exp dfw^ • Sft • w| , (2.5) 

for any functions Wa{t). Here tildes denote Fourier transforms, according to the convention 
that 

h{f) = I ^^'^'h{t)dt. 

Since the random processes na{t) and Saif) are uncorrelated, the spectral density matrix of 
the detector outputs is just the sum of those for the detector noise and for the background: 

S/.(/) = S„(/) + S3(/). (2.6) 

To derive an expression for Ss{f), two key ingredients are needed. The first is a mode 
expansion for the metric perturbation for an isotropic, stationary SB. Expressed in frequency 
space, this is 
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(s^)^.,(x,/)= /ci^n^ Y: ^An(/)e^-^^"-e^'" (2.7) 

A=+,x 

for / > 0, and (^^)/ijj(x, /) = (^^)/ijj(x, — /)* for / < 0. Here the tensors e^-" are the 
usual transverse traceless polarization tensors, normalized according to e^'"e^'" = 26 ab, 
and / (PQn denotes the integral over solid angles parameterized by the unit vector n. The 
coefficients SA,n are random processes which satisfy 



{sA,nif)sBAf'y) = SAB5\n,m) 

X 6if - n ^ n,if) (2.8) 

and 

{sAAf)sB,min)=0, (2.9) 



for /, /' > ||2T|. Here 5^(n, m) is the delta function on the unit sphere. 

The second ingredient is the expression for the response of the ath detector to the back- 
ground. This is 

s,(t)=d, :(SB)h(x„t), (2.10) 

where x^ is the position of the detector, the : denotes a double contraction, and is a 
symmetric tensor that characterizes the detector's orientation (its polarization tensor). If 
the arms of the detector are in the directions of the unit vectors 1 and m, then = 
(1 (g) I - m (g) m)/2 [||]. From Eqs. (^) and (pTToD we obtain that 



5a(/) =E / rf'^^ni^,^(n)s^,ne2-/"-^ (2.11) 

where F^{n) = da : e^'" are the detector beam pattern functions. Inserting this response 
function into an equation analogous to Eq. (^.4|), and using Eqs. (|2.8|) and (P^), we obtain 



5.(/)., = ^^^7a.(/), (2.12) 



where |^| 
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X exp [27ri/n ■ (x„ - Xf,)] . (2.13) 

The functions ^ah are the overlap reduction functions discussed in Sec. |. It can be seen 
that their value is unity for coincident, aligned detectors. Below when considering a single 
pair of detectors we shall write 7a6(/) simply as 7(/). 

III. THE OPTIMAL PROCESSING STRATEGY 

In this section we describe the optimal method for filtering the detector outputs when 
searching for the SB. We discuss the two detector case in subsection [111 A| the filtering 



method in this subsection is implicit in the formulae of Ref. [ffBl. A detailed proof that 



the method is optimal is given in Appendix A. In subsection [1IIB| we discuss the effects 
of correlated sources of noise, and explain why correlation measurements between detectors 
at widely separated sites yield much better upper bounds on ^g{f) than correlation mea- 
surements between detectors at one site. Next we describe the modifications to the filtering 
method necessitated by correlated noise, in subsection [III C. 



A. General description 

To measure the stochastic background one really needs to measure the spectral density 
matrix of the detector outputs Sh{f)- Now there is no way in principle to separate out the 
portions of Sh{f) due to detector noise S„(/), and due to the SB Ss(/). If we had only one 
detector, we could only conclude that Ss{f) < Sh{f)- From Eq. ([2.12| ) and using ■jaaif) = 1, 
this would give an upper bound on Qg of 

W)<2.5xlO-«(M|l](_i_],,|. (3,, 



Here /i„ = \J fSh{f), which is projected to be ^ 10 '^^ for LIGO, even at an advanced stage 
|l|]. The quantity hf^ is the Hubble constant scaled to the value of 75kms~^Mpc~^. 
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It is unlikely that Qg will be as large as the value in Eq. ( |3.1| ). However, if it does happen 
that Qg ^ 10^^ in the LIGO/VIRGO waveband, then the SB induced noise may dominate 
over the other sources of detector noise at some frequencies, and may ultimately constrain 
the amount of information that we can extract from burst gravitational waves. In this paper 
we shall from now on assume that Qg is small, and restrict attention to measurements made 
using two or more detectors. 

With two or more detectors, one takes advantage of the fact that the sources of noise in 
each detector will be independent. This will be the case for detectors at widely separated 
sites, because sources of noise that are correlated between the detectors on timescales of the 
order of the light travel time between them are expected to be insignificant, or if not they 



can be monitored and compensated for |17|. When correlated noise is unimportant, then 



the off-diagonal elements of S„(/) will be very small, so that 

Ss{f)ab-Sh{f)ab for a ^6. (3.2) 

By measuring these components we can gain information about the SB. One does this by 
cross correlating the two output streams |]TU|. One takes each detector output ha(t), a = 1, 2, 



and constructs using an optimizing linear filter K{t) the quantity Ha{f) = K{f)ha{f). The 
purpose of this filter is essentially to suppress the signal at those frequencies at which the 
detector noise is strong, and it is given by [cf. Eq. ([A51|) below] 

^ (3.3) 

For coincident, aligned detectors the next step is simply to integrate Hi{t) against H2{t), 



see, e.g., Ref. [|^. For noncoincident detectors, however, a different strategy is necessary. 



One first constructs the correlation with time delay r, 

Y{r) = dtHi{t + T)H2{t), (3.4) 

J-t/2 

where f is the observation time, typically of the order of a year. Then one calculates the 
weighted average 
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Y = r dTL{r)Y{T), (3.5) 

where ti is the hght travel time between the detectors, and L{t) is a weighting function 
which must be carefully chosen for each detector pair in order to maximize the sensitivity. 
Roughly speaking, this smearing of the cross-correlation compensates in some measure for 
the phase lags between the detectors which were discussed in Sec. |. The quantity Y will 
then have a truncated gaussian distribution (i.e. gaussian but restricted to positive values) 
with signal to noise ratio given by Eq. (|1.2|) ||T8| . 



In fact the sliding delay function L{t) is just the Fourier transform of the overlap re- 
duction function [see Eqs. ( A50 ) and ( |A51 )]. In Appendix B we give an analytic formula 
for L{t), and we show in Fig. ^ the sliding delay function that will need to be used for the 
LIGO pair of detectors. 

In order of magnitude, the 90% confidence upper limit that can be placed on Qg by 
cross-correlating between the detectors is |]TU 



n^-"^ ^ (3.6) 



Here fio ~ lO"*^ is the upper bound ( p.l|) obtainable from one detector, and A/ is the 
bandwidth of the measurement. [If no bandpass filtering of the data is carried out. A/ will 
be roughly the width of the peak of the function l/(/^5'„); and if filtering with a bandwidth 
A/ is used, then the domain of integration in Eq. (|1.2| ) must be suitably restricted]. 



B. Effects of correlated noise. 

Up to this point we have assumed that the intrinsic detector noise is uncorrelated between 
different detectors, i.e. that the matrix S„(/) is diagonal. We now relax this assumption and 
consider the effects of correlated noise. It is planned for LIGO/VIRGO to ultimately have 
two or three detectors per detector site, perhaps optimized for different types of gravitational 
wave sources. Thus, there will be two possible types of correlation measurements: intrasite 
measurements at one site, and intersite measurements between detectors at different sites. 
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As previously mentioned, it seems very unlikely that there will be any significant correlated 
noise in intersite measurements, and so it will only be important for intrasite measurements. 

We now consider what information can be extracted from intrasite correlations. The 
quantity we can measure is the sum S„(/) + Ss(/), the two terms of which are in principle 
indistinguishable Hence, we can draw inferences about Ss(/) only if we have some 



information about the correlated noise. The most obvious such a priori information is the 
fact that S„ is positive hermitean matrix. It follows that if we measure the total power 
spectral density matrix to be Sh, then < S^. However, this only tells us that ^Ig ^ 10~^, 
[cf. Eq. (|3.1| ) above], since the inequality applies to the total matrix, and the off-diagonal 
elements are small compared to the diagonal elements. In particular, it is not true without 
further assumptions that 

Ss{f)ab < I Sh{f)ab I for a ^ b, (3.7) 

(as is implicitly assumed in Ref. jl^]), since Sn{f)ab may be negative (or complex). A 
complex value of Sn{f)ab corresponds physically to sources of noise that excite two detectors 
with a certain preferred phase lag between them. 

In Appendix A, we derive the 90% confidence upper limit that one can place on Qg in 
a bandwidth A/ using intrasite correlations, which generalizes Eq. ( p.6| ). If we define the 
noise correlation coefficients by 

/7 Sn{f)ab /Q qN 

Y Sn{f)aa Sn{f)bb 

then the result depends on (i) the measured value C of C (determined from Sh), and (ii) the 
assumed a priori maximum value Cmax of |C|. In order of magnitude we find [cf. Eqs. 



and (|MD|) below] 

^ (3.9) 



Since there will be various unknown sources of weak correlated noise, it is not appropriate 
to choose a very small value of Cmax- Hence, the upper bound ( |3.9D will be much worse than 
the bound ( |3.6|) obtained from intersite correlations. 
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We also consider the assumption that the correlated detector noise always excite different 
detectors in phase, but can be arbitrarily large in magnitude, so that 

SnUlab >0 for a y^b. (3.10) 

In this case intrasite correlation measurements cannot be used to detect the SB, but can only 
be used to place upper bounds on its magnitude. The resulting upper limit [cf. Eq. ( [A3 81 )] 
is given by Eq. ( |3l9| ) with C^ax = 0. Thus, the intersite correlations will still give better 
bounds unless the actual amount of correlated noise in the detectors satisfies 

C<-^ (3.11) 

which is ~ 10^'^ for a year-long measurement with a bandwidth of ~ 50 Hz. Since it is 
not clear that either of the conditions ( |3.10| ) or ( |3.11| ) will be appropriate, we henceforth 



consider only intersite correlations. 

C. Filtering intersite correlation measurements. 

In Appendix A we show how to optimally filter the intersite data when there is a detector 
network with several interferometers per detector site. In this case the detector noise matrix 
S„(/) will have a block diagonal form, with each block corresponding to a site. Let 
be the subblock corresponding to the Ath site. The strategy is essentially to measure 
the off-diagonal blocks of Sh{f) (i.e. the intersite correlations), and to use these to obtain 
information about the SB. The resulting SNR squared is then given [cf. Eq. ( |A52D below] 
by a simple generalization of Eq. ( PT^ ) : 



VSvr/ Jo ' p 



too 

2f df- 
Jo 



Here the sum is over the pairs of sites, and 7ab is the overlap reduction function for any 
detector at site A together with any detector at site B. The quantities S^^^\f), which we 
call the effective site spectral noise densities for correlations between sites, are given by 
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-1 

(3.13) 

They will be real since the matrices are Hermitean. In the case where the off-diagonal 
elements of are all equal to S'^{f), and the common value of the diagonal elements is 
Snif), then we obtain that 

sf\f) = ^sM-) + {i-j^)s'M), (3.14) 

where is the number of detectors. This formula shows that the use of detectors instead 
of a single detector at one site reduces that sites effective noise density by a factor of 1/A^, 
if the detectors noise sources are independent {S'^ ~ 0). If the noise sources are strongly 
correlated, however, so that S'^ ~ Sn, then there is no significant reduction in the effective 
noise. 

IV. THE OVERLAP REDUCTION FUNCTION 

In the preceding sections we have seen that for a pair of interferometric detectors, the 
overlap reduction function 7(/) characterizes the dependence on the detector separation and 
orientations of both the correlation matrix ( p.l2|) , and the broadband sensitivity to the SB 
( |1.2|) . We now give an analytic formula for this function and discuss some of its properties. 

A. The general formula for 7(/) 

We first define the variables which describe the orientations and separation of a pair of 
detectors. Let L be the line joining the two detectors and Pi (P2) be the plane formed by 
the arms of the first (second) detector; see Fig. ^ The variables we will use are: (i) The 
distance d between the detectors; (ii) The acute angle Pi between L and Pi; (iii) The angle 
(Ji between the projection of L onto Pi and the bisector of the two arms of the first detector; 
(iv) Corresponding angles /?2 and a2', (v) The angle x between L and and the intersection 
of Pi and P2. The directions (clockwise or anticlockwise) in which ai and (T2 are chosen 
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to be positive are unimportant, as long as the conventions for ai and for a2 coincide as 
A;/52 0. Let 6 = (o"i — o'2)/2 and A = (cti + (72) /2. In Sec. | we have called these the 
relative rotation and total rotation angles respectively, since for (3i = (32 = 0, the angle 6 
is half of the relative rotation of the detectors while A measures the average rotation with 
respect to the line joining them 



In general the overlap reduction function as given in Eq. (|2.13|) will depend on all of the 
variables P2, A, x, and on the phase lag a = 27ifd between the detectors, where / is 
the frequency. Now for terrestrial detectors x = '^/'^ ^ind (3i = (32 {= P say), which is also 
the angle subtended at the center of the earth by the detector pair. For this case we derive 
in Appendix B the following formula: 

7(/) = cos(45)ei(«,/5) + cos(4A)e2(«,/3), (4.1) 

where the functions 61 and G2 are 

ei{a,(3) = cos\(3/2)g,{a) (4.2) 

and 

e2{a,(3) = cos\(3/2)g2{a) + g3{a) 

-sm\p/2)[g2{a)+giia)]. (4.3) 



The functions gj{a) are given in Eqs. (pil|) , ( |B12| ) and (pi3| ) of the Appendix, and are all 



linear combinations of the functions sin(a)/Q;" and cos(Q;)/a" where 1 < ri < 5. 

Several properties of the overlap reduction function 7(/), which were discovered by Chris- 
tensen and others from numerical calculations |T^ , p!6|jr?| , can be read off the formula (^H]). 
First, there will always be frequencies / for which 7(/) vanishes, and correspondingly near 
which the narrow-band sensitivity of the detector pair to the SB is very poor. For de- 
tectors that are less than a few thousand kilometers apart, the first null frequency is at 
/i ~ (70IIz)(3000km/(i), irrespective of the detector orientations (see below). Second, the 
reduction function falls off like 1/f when a ^ 1, or equivalently when f ^ fi- Hence the 



90% confidence limit that we can place on ^g{f) scales like l/{fd) for large d. However for 
detectors which are ~ 3000 km apart, the phase lag a is of order unity for typical detector 
frequencies, and hence the variation of the sensitivity with distance is more complex than 
simply scaling like 1/d (see Sec. |VT| below and also Fig. |I] above). 



B. Some special cases 

To understand the behavior of 7(/) as a function of 6, A and /5, it is useful to consider 
some limiting cases. For coplanar, coincident detectors, d = (3 = 0, and using gi{0) = 1 
and 5'2(0) = (73(0) = we find that 7 = cos(45). This is just what we would expect 
physically. Consider, for instance, detectors which are rotated with respect to one another by 
45°. Vertically incident gravitational waves will couple to the two detectors via polarization 
components that are orthogonal and hence statistically independent; thus these waves give no 
contribution to the cross-correlation. Waves which are incident from non-vertical directions 
do give rise to correlations between the detectors, but when we average over all incident 



directions the total correlation vanishes |26|. It is also clear from rotational symmetry that 
there should be no dependence on the total rotation angle A. 

Now suppose that the detectors are still coincident but no longer coplanar, so that 
/9 7^ 0. Equivalently, suppose the detectors are separated by a distance d but evaluate 7(/) 
at frequencies / <^ 1/d. The result is 

= cos^(/3/2) cos(45) - sm\(3/2) cos(4A). (4.4) 



This equation agrees with Eq. (3.9) of Ref. after making an appropriate change of 
variables. As we would expect, the dependence on A is small when f3 is small; it is smaller 
than the dependence on the relative rotation angle 5 by a factor of tan^(/3/2) (~ 3 x 10~^ 
for the planned LIGO detectors; see below). 

A third limiting case is when /3 = but d^ 1//, so that the detectors are coplanar but 
effectively separated. In this case. 
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7(/) = cos(4(5)(7i(«) + cos(4A) {52(a) + g-sia)} . (4.5) 

Plots of the functions gi{a) and (72(a) + (73 (a) are shown in Fig. ^. We see that for a < 4, 

the cos(45) term dominates, just as above for coincident detectors for small (3. When the 

phase lag a is large, the functions in Eq. ( |4.5|) can be approximated as go{a) ~ 5 sin a/ (16a) 

and 5'2(a) + (73 (a) ~ —(70(a), so we obtain 

5 sin a , , , , , ^ , 

7 / ~ 5 sm 2ai sm 2^2 . 4.6 

What is happening physically here is that the cross- correlation is dominated by modes whose 
propagation vectors are nearly parallel to the line joining the detectors. This can be seen 
by applying a stationary phase argument to the integral ( p.l3| ). Hence it is reasonable that 



the overlap reduction function should have separate factors for each detector that depend 
on how they are oriented with respect to the line joining the detectors. In this case the 
dependence of the overlap reduction function 7 on the total rotation angle A is strong. 

In summary, the dependence on the relative rotation angle 6 will dominate over the 
dependence on the total rotation angle A unless two conditions are satisfied: (i) /3 ^ 7r/4, 
and (ii) a > 1, i.e. d > c//„, where /„ ~ 100 Hz is a typical frequency at which the detector 
is sensitive. In reality of course /? and d are not independent; they are related by 

d = 2sm{P/2)rE, (4.7) 

where ve is the radius of the earth. Due to the coincidence that the quantity fnfE/c is of 
order unity, conditions (i) and (ii) above become satisfied at approximately the same value 
of /3. 

For a pair of terrestrial detectors, as (3 and d are increased from zero, the following three 
effects occur: (i) 7(/) decreases because d is increasing, as discussed above (ii) 7(/) decreases 
because (3 is increasing. From Eq. ( [4.4| ) we see that the maximum over all orientations of 
the value of 7(0) is 1 — siv? (3/2. Hence the effect of increasing [3 contributes typically a 
factor of at most 1/2 to the decrease in 7(/). (iii) The dependence of 7(/) on A increases 
to become as strong as the dependence on 5. 
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V. OPTIMIZATION OF DETECTOR ORIENTATIONS 



We now turn to the issue of how to optimize the orientations of a single pair of detectors 
so as to detect or measure the SB. We assume that the detector's noise curves are the same, 
so that Sn{f)aa = Sn{f)bb = Wc cousider separately the cases of maximizing the 

narrowband sensitivity to the SB in the vicinity of some frequency /, and of maximizing the 
overall broadband sensitivity. 

Now if we rotate a detector through 90°, then its polarization tensor d [cf. Eq. ( ^.1U| )] 
will have its sign flipped. Hence |7(/)| is a periodic function of both of the angles (Xi, Ui 
(or 6, A ) with period 90°. In the following discussion, all values of these angles and all 
equations involving them should be treated modulo 90°. 



A. Narrowband sensitivity 



From Eq. ( |1.2|) we see that the 90% confidence hmit that we can place on ^g{f) in a 



small interval of frequency of width A/ centered about /, using a measurement of duration 
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.r'(/.A/).(^)/^^. (5.1) 

To minimize this we need to maximize the quantity |7(/)|. From Eq. ( [4.1|) it follows that 
the maximum value over all orientations of |7(/)| is |Gi(/)| + |02(/)|- 

Now if the signs of 0i and O2 are different, then this maximum will be achieved at 
0-1 = (T2 = 45° (i.e. either 6 = 45°, A = 0° or 5 = 0°, A = 45°). This corresponds to 
each detector having one arm parallel to the line which joins them. We will call this type 
of orientation configuration I. This will be the optimal orientation at zero frequency, since 
0i(/ = 0) = cos^(/?/2) > 0, and e2(/ = 0) = -sin^(/3/2) < 0. It will also be optimal 
at high frequencies {fd ^ 1), since from from Eqs. ( [4.1|) and ( pil|) - ( |B14| ) we see that 
0i(«) ^ -02(a) for a > 1. 
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If, on the other hand, the signs of Bi and 02 are the same, then we see from Eq. ( |4.1|) 
that |7(/)| will be maximized at o"i = o"2 = 0° (i.e. either 6 = A = 0°ot6 = A = 45°). In 
this case, which we will call configuration II, each detector arm makes an angle of 45° (mod 
90°) to the line joining the detectors. This will be optimal, e.g., for 10 Hz ^ / ^ 40 Hz when 
p = 29° (the LIGO value; see below). 

Now suppose that the detector pair is in either configuration I or configuration II. Then 
for some frequencies the orientations will be optimal, while for some other frequencies they 
will not be. It is useful to calculate, for a given frequency /, the factor 7l{f) by which |7(/)| 
is reduced if the wrong configuration out of I and II is chosen. Note that choosing between 
I and II is equivalent to fixing the relative rotation angle 6 to be zero, and then choosing A 
to be either 0° or 45°. Hence TZ measures the sensitivity of 7(/) to the total rotation angle 
A when the detectors are parallel. From Eq. ( [4.1| ), we obtain 



7^ = 



Il01l 




|02| 




|0i 


1 + 


|02| 





(5.2) 



At low frequencies {a = 27Tfd <^ 1) this is approximately 

K(/-0)«(l + ^)". (5,3) 

where u = sin(/3/2). Thus for small values of (3 we have TZ{f ^ 0) ~ 1, and whether 
configuration I or II is chosen is relatively unimportant. 

At high frequencies (a ^ 1), however, we find that 7^(/) asymptotes to a constant T^oo 
independent of frequency, where 

which can be quite small. [T^oo ~ 0.2 for the LIGO detectors; see below.] Hence the choice of 
A will have a large effect on the high frequency sensitivity to the SB, especially for detectors 
which are close together for which u is small. 
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B. Broadband sensitivity 



It is quite likely that the SB will be so weak that it will only be detectable (if at all) 
by integrating over a broad range of frequencies, e.g. from ~ 20 Hz to ~ 70 Hz. Hence it 
is more important to optimize the overall broadband response ( |1.2| ) of the pair of detectors 
than to optimize the narrow band sensitivity at some frequency. We now determine how to 
do this if we assume that the spectrum is approximately constant over the relevant 

frequency range - the so-called "Zel'dovich spectrum" . 

By inserting Eq. ( [4.1|) into Eq. ( p..2|) , we find that the SNR squared after optimal filtering 



IS 



Q2 

— = A{p) cos^(4(5) + 2B{(3) cos(45) cos(4A) 



+ C(/3)cos2(4A), (5.5) 



where 



and A and C are given by similar formulae with Oi02 replaced by 6f and 02 respectively. 
Note that A, B and C depend on (3 in two distinct ways, since from Eq. ( [4.7| ), 

ei,2(/) = ei,2(« = 47r/rssin(/?/2),/5). (5.7) 



The explicit formula (|5.5| ) makes it easy to determine how to optimally orient the detec- 
tors at a fixed value of /3, and also how the SNR at the optimal orientation varies with (3. 
The SNR ( ^.5|) is maximized at | cos(45)| = | cos(2A)| = 1, and the optimum orientations 



are just configurations I and II again, applying when 5 < and i? > respectively. 



VI. IMPLICATIONS FOR LIGO/VIRGO 



The planned LIGO detectors in Hanford, Washington and Livingston, Louisiana will 
have (3 = 27.2°, 6 = 44.9° and A = 28.2° (the last two mod 90°) [gT]. A graph of 7(/) for 
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this configuration is shown in Fig. |^. In Fig. ^ we show what the functions Gi and G2 look 
like at the value (3 = 27.2° appropriate for LIGO. For any possible LIGO orientations, 7(/) 
will be a linear combination of these two functions as in Eq. ( [4.1| ). From Fig. |^ we can see 
that there will be a zero of 7(/) at / ~ 70 Hz, whose position is relatively insensitive to the 
orientations chosen. 

The noise power spectral density in the LIGO detectors after some years of operation 
might be roughly 

Snif) = max [Sr,^{f/ fm)~^,Smif/ fmY\ , (6.1) 

where Sm = 10~^^ Hz~^ and fm = 70 Hz. For frequencies less than ~ 10 Hz the noise will 
be effectively infinite. Eq. (|6.1|) is a crude analytic fit to the noise curve estimate for the 
"advanced detectors" given by the LIGO team in Ref. [^. A more detailed model of the 
noise is unnecessary because of the uncertainty as to what the actual noise levels might be. 
We now estimate what the effective site noise level S^^^^ [cf. Eq. {^J§] might be if - 



after an upgrade of the LIGO facility - there are eventually two or three interferometers at 
each site, as is planned. The noise spectral density S'„(/) in each detector will be a sum of 
the form ]l],|ig 

including contributions due to seismic noise, thermal noise, photon shot noise and residual 
gas noise, among others. By contrast, the typical value 5*^ of the off-diagonal elements of 
the spectral density matrix will be approximately a sum of terms that are due to sources of 
noise which are strongly correlated between different interferometers at the same site: 

^ 'S'seismic ~l~ 'S'gas- (6-3) 

Since S'seismic + S'gas ^ Sn whcu / ^ 10 Hz for the advanced detectors of Ref. [|1|, we find 
from Eq. ( ^.ltj| ) that to a good approximation for LIGO, 

Sf\f)-'^j:Sn{f)-a'. (6.4) 
a 
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Here the sum is over the different detectors, and a = 1 corresponds to the primary, broadband 



interferometer with noise curve ( |6.1| ). Now the second and/or third detectors at each site 
will most likely be specialized ones such as dual-recycled detectors |^, which have high 
sensitivity only in some narrow frequency band. Thus in the relevant frequency band of 
20Hz < / < 70 Hz [cf. Fig. g below], it is likely that S^^^\f) ^ S'„(/)n. If, instead, 
three identical broadband detectors are operated at each site (which is unlikely initially), 
then Si'^\f) ^ S„(/)ii/3, and the 90% confidence upper limit estimated below should be 
divided by 3. 

If we now make the plausible assumption that the SB has approximately constant ^g{f) 
over the above waveband, then by inserting Eq. ( |6.1|) into Eq. ( |3.12|) we can calculate the 
90% confidence upper limit that we can place on Qg. Using the planned LIGO positions and 
orientations we thus obtain 



(B)f^(90%) ^ i_g5 ^ 



POO 

2f df 
Jo 



7(/) 



2 



' - \ -1/2 



-1/2 



= 5.1 X 10-'' h^,'Ni^-^^j , (6.5) 

where h^^ is the Hubble constant scaled to the value 75kms^^Mpc^^, f is the observation 
time, the prefix (B) indicates broadband, and = 5^/10^^^ Hz~^. Note that this is the 



upper bound obtained from intersite cross-correlations, since as argued in Sec. [Ill B| , it will 
not help much to include the information from intrasite correlations. 

This bound ( |6.5| ) is a reasonably conservative estimate of what the LIGO sensitivity 
might be when the advanced detectors are operating. However the ultimate noise levels are 
largely unknown. This is because the experimental techniques and technology, which are 
presently somewhat far from the advanced level, will be steadily improving. To obtain an 
upper limit to we take S^^^^f) ~ Sgssif), which is a source of correlated noise that 

cannot be changed once the beam tube is built. Using S'gas ~ 2.2 x 10~^°Hz~^ |28| gives as 
a rather firm upper limit for the LIGO sensitivity 

(B)^(90%) > 2 X 10-i=^iV/i7t (y^) • (6.6) 
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We now consider how would change if we were to vary the angles /?, S and 

A. The integral (|5.6|) is positive for (3 = 27.2° with the noise power spectral density ([Ol). 
Hence the optimal configuration of the detectors is configuration II, i.e. at 5 = A = 0° 
or 45°. However LIGO already has 6 ~ 45° (corresponding to parallel detectors because 
cTi — = 2(5 ~ 90°). Also, as we have discussed in Sec. the dependence of 7(/) on A 
is very weak for / < 100 Hz and at /5 = 27°, see Eq. ( [4.1| ) and Fig. ^. Because of these two 
facts the expected LIGO sensitivity at its planned configuration of 6 = 44.9° and A = 28.2° 
is only ~ 3% less than the optimal sensitivity which would be attained at 6 = A = 45°. 
This follows from Eq. where we find that B/A = 0.045 and C/A = 7.8 x 10^1 

The European VIRGO and GEO detectors, however, are currently planned to have 2S = 
(Ti — (72 ~ 45°. Such a configuration will allow VIRGO/GEO to extract information from 
both polarization components of burst and periodic gravitational waves, but will severely 
reduce the sensitivity to the SB. For these detectors /5 = 9.75°, and we find from Eq. ( ^.6|) 
that B/A = 4.0 x 10"'' and C/A = 1.7 x 10~^. In Fig. we plot the broadband signal to 
noise ratio ( |5.5| ) as a function of 6 (assuming that A is optimally chosen); it can be seen 
that at 6 = 45°, the sensitivity to the SB is reduced by about three orders of magnitude. 
An intermediate choice of 5 ~ 40° would give a broadband sensitivity to the SB comparable 
to that of LIGO, while still allowing the detectors to access two polarization components of 
the gravitational- wave field that are largely independent. 

The total rotation angle A is unimportant for the broadband sensitivity to the SB in 
the case of detectors which are sufficiently close to each other [but not unimportant for the 
narrow-band sensitivity, cf. Eq. (|5.4| ) above]. However, for detectors which are on separate 
continents, A will become important. For example in Fig. ^ we show the functions Oi(/) 
and 02(/) for P = 79.5°, the value appropriate for one of the two LIGO/VIRGO cross 
correlations. Since 62 is typically large compared to Oi in this case, varying A will have a 
large effect on the sensitivity. This is confirmed by Eq. {^^), where we find B/A = —0.22 
and C/A = 18.3. 
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Next we consider how the broadband sensitivity to the SB varies with the angle (3, 
assuming that the orientations are optimally chosen. Recall [cf. Sec. |V[] that this means 
either configuration I or II is chosen. In Fig. |I] we plot the SNR ( |5.5| ) as a function of (3, for 
both configurations I and II, again assuming the noise curve ( |6.1| ). It can be seen that the 
sensitivity falls off rapidly when the detectors become more that a few thousand kilometers 
apart. The sensitivity for the LIGO separation is already a factor of ~ 4 worse than that 
for coincident detectors. Hence if there is ultimately a worldwide network of detectors 
in America, Europe, Japan and Australia, only the cross-correlations between proximate 
detectors like LIGO/LIGO or VIRGO/GEO will be important. 

The various conclusions that we have drawn are relatively insensitive to the the detailed 
properties of the detector noise curve (|6.1| ) that we have assumed (except for the value of 
the minimum frequency). This is because the integral ( |1.2|) is dominated by contributions 
in the narrow frequency band 20 Hz ^ / ^ 70 Hz; see Fig. ^ At these frequencies thermal 
noise dominates over photon shot noise [Q, so using sophisticated optical configurations of 
the interferometer such as dual recycling PD| to reduce the shot noise will not help much. 
It thus will be important for the detectors to have good low frequency sensitivity for the 
purposes of placing upper limits on the SB. 



VII. CONCLUSION 

To place upper bounds on the strength of the gravitational stochastic background, or 
perhaps to detect it, is but one of the aims of LIGO and other detector systems. Moreover, 
this background may well be very weak compared to waves from astrophysical sources; 
consider by analogy the weakness of the electromagnetic stochastic background in the visible 
region of the spectrum. However, these waves if detected would be amongst the most 
interesting that the detectors would see. 

The sensitivity to these waves is therefore just one of various factors that need to be taken 
into account when choosing the fixed orientations of interferometric detectors. Nevertheless, 
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in order to be able to make a wise choice, it is important to know the effect for the sensitivity 
of choosing this or that orientation. One of the main results of this paper is the simple 
expression ( |B5[ ) describing this dependence. In particular, for frequencies such that the 
phase lag between the detectors is ^ 1, the sensitivity is determined solely by the overlap 
of the polarization tensors of the two detectors. At higher frequencies the direction of the 
vector which joins the detectors also becomes important. 

We also determine the orientation of a pair of interferometers that will optimize the 
sensitivity of the pair to the SB, and show that the orientations which have been chosen for 
the LIGO detectors are close to optimal. 
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APPENDIX A: STATISTICAL FOUNDATIONS 



1. Overview 



In this Appendix we justify the method of optimal filtering described in Sec. |Tl| and the 



resulting signal to noise ratio ( 3.12 ). We do this by deriving the probability distribution for 
the spectrum Qg given the detectors' output h{t), p[Qg\h{t)]. 

First we transform from a continuous to a discrete description of the measured data. In 
the body of this paper we have treated the output of the detector network as a continuous 
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vector random process h(t). However, a real, discrete measurement will be of a finite du- 
ration f, and will contain frequencies only up to some maximum frequency determined by 
the time resolution of the sampling. In other words the output of the detectors will consist 
of the numbers haj = ha{tj = tstart + J A) for 1 < a < Ud and < j < N, where rid is the 
number of detectors, tstart is the starting time, and A is the time resolution, of the order of 
10~^s. The number of samples per interferometer is = f/A. We denote by X the vector 
of numbers haj, which we assume to have a multivariate gaussian distribution with some 
variance-covariance matrix S'. This matrix is essentially the correlation matrix (|2.2| ): 

K^bJ=Ckit^-t,U. (Al) 

However if we take a finite Fourier transform, which amounts to making a change of basis 
in the space of vectors X, then S' corresponds instead to the spectral density matrix (p7^). 

Now from Eq. ( |2.6| ), we have S' = + S(,(f2), where the contribution from the SB 
depends on the spectrum ^g{f ) as in Eq. ( p.l2| ). [Because of the finite frequency resolution 
of order ~ 1/f, we represent the function ^g{f) as a finite vector fl.] What we want to do 
is to extract information about f2 from a measurement of X. 

There are two general approaches to the task of quantifying the information obtained, 
in situations of this sort. First, one can in principle compute the probability distribution 
function (pdf) for f2 given the measurement X, p(f2|X). This pdf then contains complete in- 
formation about the measurement. However, the calculation of the pdf p(ri|X) is frequently 
difficult in practice; and so one has to resort to instead calculating estimators (statistics) 
f2(X) which are functions of X, chosen so that the pdf for their values given some value of 
fl, p{fl\fl), is peaked near f2 = f2. There are standard criteria for choosing such estimators, 
see for example Refs. pO| , |3T |. 



Now suppose that, instead of one measurement of X, we have n measurements 
Xi, . . . , X„. A standard statistical result is that in the Cramer- Frechet-Rao limit of n ^ oo 
PT| , P2[| , the two approaches discussed above yield the same unique result. More precisely, 
the pdf p(r2|Xj) becomes a gaussian centered at f2mi(Xj), where fl^\ is the so-called max- 
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imum likelihood estimator of fi. The variance-covariance matrix of this probability 
distribution depends on the only through Cl^\(Kj), Eq = SQ(f2jni). Conversely, the 
pdf p(fimi|f2) of the estimator fimi given some value of f2, becomes a gaussian centered at 
f2 with width 5]n(r2). Thus the two probability distributions, which are conceptually very 
different objects, become effectively the same, and all one really needs to calculate is the 
variance-covariance matrix Sn(ri). 

However, this simplifying Cramer-Frechet-Rao limit does not apply in a straightforward 
manner to the calculation of p{fl\^) for measurements of the SB; there are a number of 
subtleties and differences from the usual situations discussed in Refs. 0,^. First, we have 



only one measurement of X instead of a large number n of measurements. Nevertheless, 
something like a "large number of measurements limit" does apply, which we discuss farther 
below. Second, one needs to address the issue of distinguishing between the contributions 
to S' from the detector noise and from the SB, since both are unknown a priori. Third, 
the CFR limit only applies to the extent that our a priori knowledge of the variables being 
measured is unimportant. However, our a priori knowledge about the correlated detector 
noise in intrasite correlations is important in the calculation of p(ri|X), and consequently 
this pdf is not approximately gaussian [cf. Eqs. ( [A21| ) and ( [A3 3D below]. 

To resolve these complications, we now derive from first principles an approximate ex- 
pression for p(f2|X). We also identify the conditions under which the approximations we 
make are valid, and show that they will all be satisfied in the LIGO/VIRGO context. We 
assume throughout that the effect of the SB on the detectors is small compared to the 
detectors intrinsic noise, so that from Eq. (|0|), ^g{f) ^ 10~^. 

2. Calculation of p(n I X) 

We start by considering the frequency resolution of the measurement of Qg{f). In a finite 
measurement of length f there are roughly N = t/A independent frequencies fj, and we get 
what amounts to one measurement of each Qj = Qg{fj). Now if f is doubled, we double the 
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number of measurements, but we also double the number of variables measured. Clearly, 
the Cramer-Frechet-Rao limit of repeated measurements of the same variables cannot be 
attained in this context without further assumptions. Below we shall argue that in the 
LIGO/VIRGO context, to a good approximation, 

Ch{T)^0 for r>r„ (A2) 

for some correlation time ^ f [Q. In this case the unknown quantities to be measured 
are C/i(jA) for < j < = t^/A, or equivalently Sh{f) at n different frequencies; and the 
number of measurements of each of these variables is ~ N/n = t/tc^ 1. 

In order to explain Eq. (S3), we fix a value of Tc, and decompose C/i(r) into 



C,,(r) = C;,(r)e(r, - |r|) + C,(r)e(|r| - rj, (A3) 

where 9 is the step function. Let Sh = + be the corresponding decomposition in 
frequency space, so that S^(/) is Sh{f) averaged over frequency scales of the order of l/r^. 
Then from Eq. ( |2.6| ) there will be two contributions to S^(/): 

S^(/) = S>(/) + S>(/). (A4) 

The first term here will be small if Qg{f) is smooth over frequency scales ~ I/tc [cf. Eq. ( |2.12|) 
], which will be the case for currently conceived of SB sources when = 100 s, for example. 
The second term in Eq. (|A4|) will also be small for this value of Tc, except near isolated 
frequencies corresponding to high-Q resonances [Q, the effect of which can be neglected. 

We now calculate the pdf for Q using Bayes rule; by virtue of the condition ( [A2D , 
something like the CFR limit will apply. Let S' = X (g) X, which is the maximum likelihood 
estimator of S'. This is a Nud x Nrid matrix, where nd is the number of detectors. Similarly 
the matrices and S'^ are of dimension Nud, and are given by formulae analogous to ( [All) . 
However because of the condition (|A2|) the independent variables in S' can be formed into 
a smaller matrix S, which is given by 

Eaibj =Ch[{t-j)AU (A5) 
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for < i, j < n = Tc/A. Note that since C/i(— r) = C/i(r)'^, the matrix S contains only 
n'^n independent variables. We define the matrices S„, and S analogously in terms of 
C„(r), Cs(r) and the estimator 

\ ftsta.it+T 

ChiT)ab = - dt K{t + T)h{t), (A6) 

^tstart 

which is defined for |r| < r^. The estimator C/i is the quantity that will be measured; it 
is a sufficient statistic for O, and moreover its discrete counterpart S(X) is the maximum 
likelihood estimator of S. 

The joint pdf for S„ and given X is 

1 



2 ^ ' 



(A7) 



p(S„,S,|X) (X pi°)(S„)pio)(s,)exp 
where 

A'(S') = lndetS' + S' : (A8) 

and p^^^ and p^'''' are the pdfs that represent our a priori knowledge. From the distribution 
(|A7|) we obtain the pdf p(0|X) for f2 given X by integrating over S„: 



p{n\X)=U j rf5]„p(S„,S,(l^)|X), (A9) 

where A/" is a normalization constant. 

The above formulae involve primed, Nud x Nud matrices. We now express them in terms 
of the corresponding, unprimed, nud x nud matrices, by using the fact that all the matrices 
will be approximately diagonal in frequency space. That is, they will be approximately 
diagonal in the indices i,j after a finite Fourier transform change of basis, but not in the 
indices a, h. For example, on a suitable basis, 

{^n)ai,j^5ljSn{fl)a,, (AlO) 



where fj = I/Tc for 1 < I < n. It is straightforward to verify using Eqs. ( |A^ ) and ( |A6D that 
35| 
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A'(S') ^ fcA(S), 



(All) 



where k = t/tc = N/n is the effective number of measurements, and 



A(S) =lndetS + S : S'^ 



For example, the second term in A'(S') can be written as 



(A12) 



1/ . v-i'-l 



S XaXjS 



2 1 d/h(/)t.S,(/)-i-h(/), 

J —CO 

r -^1 

f dfti s,(/)-i ■ S,(/) . 



(A13) 

J — oo 

Here we have used Eq. ( p.5| ), switched into the time domain, used Eq. ( |X6| ) and switched back 
| 35| . The result is just kfl : and the second term in Eq. ( |X8| ) transforms analogously. 

Now Eqs. ([A7| ), ( |A9|) , ( |A11|) and ( |A12| ) together yield a formal expression for p{Q,\X). To 
proceed further we need to invoke some approximations. The first of these is the quadratic 
approximation: we expand the likelihood function ( |A12| ) to give 



A(S) = Indet E + nrirf 
1 



+ -tr 
2 



(AM) 



where 5S = S — S = S„ + Ss(ri) — S. Clearly this approximation will be good only for 
certain values of I]„ and of fi; we discuss below the implications of this restriction. First 
we derive the conditions under which most of the probability of the pdf 

k 



p(S|X) oc exp 



:A(S) 



(A15) 



will be concentrated in the region Q where the quadratic approximation (|A14|) is valid, so 
that the normalization of the pdf ( |A15|) can be calculated correctly using Eq. ([A14| ). If we 
let Ai, . . . , A„ be the eigenvalues of S ■ S^^ — 1, then up to an additive constant, 

1 



A(S) 



E 



ln(l + A,) + 



1 + A, 



(A16) 
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From this formula it can be seen that the quadratic regime Q is given by max^ \Xj\ < e 
where e is some small number. If we demand that the total probability in Q be 1 — 5 for 
some 5 <^ 1, and assume that the a priori distribution p^^^ is not sharply peaked, then we 
find the condition 

AT > 4(logn+ |log(5|). (A17) 

This condition will be just barely satisfied for, e.g, e = 5 = 0.01, f = 10^ s, Tc = 100 s, and 
A = 10~^ s, which are values that are appropriate for LIGO/VIRGO. Now if ^ 10~^, then 
the intgeral over 5]„ in Eq. (|A9|) will be dominated by contributions from the region Q, and 
so the approximation is valid. For ft ^ 10~^, the exact probability p(0 | X) is very small, 
and so the quadratic approximation will still give a qualitatively correct result. In particular, 
the normalization of the pdf p{fl \ X) resulting from Eq. ( |A14| ) will be approximately correct. 

The next approximation involves consideration of the relative magnitudes of various 
components of the measured autocorrelation matrix I] [cf. Eq. ( |7^ ) above]. We introduce 
the following notation: for any matrix A, A" denotes the matrix consisting of the diagonal 
subblocks (in the indices a, b) of A corresponding to intrasite correlations, and A-*- = A — A" 
consists of the off-diagonal subblocks. We also define A"^ to be the matrix of diagonal 
elements (in the indices a, b) of A and A*^ = A" — A-^. Thus, the estimator S decomposes 
into + S*^ + S^, where S"^ contains the detector noises Ch{T)aa, contains the 
measured intrasite correlations, and S"*" the measured intersite correlations. Now we expect 
the contribution of correlated detector noise to to be very small, cf. the discussion in 
Sec. pi C| above. Hence from Eq. ( p.l|) , S-*- ~ eS^, where 

For example, if the SB is just barely detectable, then e ~ 10~"^. We similarly define in order 
of magnitude the small parameter e by 

S^~£S^, (A19) 
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which we expect to be of the order of 5'gas(/)/5'„(/) ~ 10 ^ or smaller if Qg < 10 ^; see 
Eq. (U. 

In the expression ( |A14| ) for the likelihood function, to leading order in e we can replace 
the factors of S""*^ by (S")"^, so that the cross terms between 61]^ and ^S" vanish. Note 
that to this order, (S"^)" = (Sll)^^ Hence, the pdf splits into a product of two factors 
which incorporate the measured intersite and the intrasite correlations. The a priori pdf for 
the detector noise in Eq. ( |A7|) can be written as 



(A20) 



since we expect S;J; ^ 0. Using Eqs. 
k r 



(IMTD , ( \KT^ ) and ([Mod we obtain 



X exp tr 



4 



(sii)-i-(s,(n)^-s^ 



(A21) 



where A/i is a normalization constant. Here 

pc{n)^ Jd^l pi°j|(Si) 

k 



X exp < — — 



(A22) 



is a function representing the information from intrasite correlations, and 5S" = S|[ + 

s,(n)ii-sii. 

Now if £ ^ e as estimated above, then most of the information we obtain from the 
SB will come from the intersite correlation measurements. In this case the factor pc{^) 
in Eq. ( |A21|) can be approximated to be constant [see Eq. ( |A33| ) below]. We discuss the 
implications of the resulting pdf for data-processing in subsection [A4| . However it may turn 
out that e ^ e, which from Eqs. ( [A IS] ) and ( |A19| ) will be the case if Qg ^ 10^^. Even if 
Qg < 10"^ it may happen that the correlation coefficient between different detectors for 
residual gas noise in the beam tube will be small compared to unity, so that in the notation 
of Eq. (|6.3D, S'^ ^ Sgas and again e ^ e. Alternatively, it may be possible to detect and 
take account of bursts of outgassing from the beam tube walls. More detailed discussions 
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of the possible magnitude of intrasite correlated noise can be found in Ref. [|T7|. Because of 
the possibility that e < e, we now derive an approximate expression for pc{^)- 



3. Information from intrasite correlations 

To evaluate the expression (|A22|) we first make a change of variables. Let C be the matrix 
of correlation coefficients 

and similarly define matrices N and G by replacing S*^ in Eq. ( [A23| ) by and 5]s(ri)ll 
respectively. Also let 

C = (S^)-^/^■S^(5]^)-V^. (A24) 
We assume that the volume element in Eq. ( |A22| ) can be written as 

rfSjl OC Pnoise(S^)pcorr(C)dNdC. (A25) 

To leading order in e the cross term between 611^ and 5S*^ in the argument of the exponential 
in Eq. (|A22|) vanishes, and so it becomes proportional to 



tr 



(n + g^-i)' + (ni/2 . C ■ + - C 



(A26) 



Using the relations C~e<^l, G~ Qg/10 ^ ^ 1, and \C\ ^ 1, and assuming that 

the pdf pnoisei^n) is slowly varying, we can approximately carry out the integral over N. 
From Eqs. (|A22D , ([K25| ) and ([X26D , the result is 



Pci^) ~ A/'2 J dC PcovviC) 



xexptrj-^ [C + G^-C]^|, (A27) 

where M2 is a constant. This expression is actually only a good approximation when \C\ <^ 1, 
but the integral ( [A27] ) is dominated by values of C close to C — G*^ which is <^ 1. 
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Now the behavior of the function ( |A27|) depends strongly on our a priori information 
about the correlation matrix C. As in Eq. (|A10|) , this matrix will be approximately diagonal 
on a frequency basis: 



where from Eq. ( |A24j ) 



abl 



-al,bj ~ SljCabI 



Sn{fl)ab 



(A28) 



(A29) 



\J Sn{fl)aa Sn{fl)bb 

Since Sh{f) is positive definite, we have \Cabi\ < 1- Moreover, the variables Cabi will be real 
whenever 



Mt + t) ribit)) = {na{t) nb{t + r)) 



(A30) 



for all t and r, and thus in particular they will be real if there are no sources of noise that 
preferentially excite one detector later than another one. If we assume that such sources of 
noise are negligible, then by inserting into Eq. ( [A27| ) the pdf 



p,orr(C) = 5(ImC)p;„„(ReC), 



(A31) 



we obtain an equation again of the form ( |A27| ), but where C and C are now understood to 
be real. [A factor of exp — tr (ImC)^/4 is absorbed into A/2.] 

Next we assume that our a priori information about the correlation coefficients Cabi is 
of the form C^in < Cabi < Cmax, so that we can take p6 | 



Pcorr{c)dc = n n 

7=1 a<b 

X 0(Cmax — Cabi) Q{CabI — Cmin) dCabI , 



where is the step function. Using Eqs. ( |A27| ) and (|A32|) and the replacement tr 
[cf. Eq. (IAT3D ], gives 



(A32) 



pc{n) oc n n 

1=1 a<b 



erf (C^i^^Vk - I3abi) - erf (CminV^ - (3abi 



(A33) 
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where 



abl 



(A34) 



If we define Sh{f) to be twice the Fourier transform of the estimator ( [A6| ) as in Eq. 
and put 

5nf 



n, 



abl 



Snifi) 



o 

abl 



(A35) 



then we find from Eq. ( p.l|) that 
4pc 



abl 



5^P\SHifl)aaShifl)bb 



10-6 



^abl — ^l) , 



(A36) 



where Qj = Qg{fi). 

We now discuss the content of Eqs. (|A33 ) - ( A36 ). If we make no assumption about 
the amount of correlated noise present, and so take Cram = — 1 and C^ax = 1, then we see 
using k ~ 10^ that pc{^) is roughly constant for < Qj 10^^. Thus we obtain little 
information. More information about the correlated noise needs to be input in order to 
constrain the strength of the SB. One assumption which may be valid is that 



Sn{f)ab > 0, 



(A37) 



which can be enforced by setting Cmm = in Eq. ( |A32|) . Equation ( |A37|) will hold if there 
are no sources of noise h(t) which contribute an amount +h{t) to the output of one detector, 
and an amount —n{t) to the output of another. As example of such a source of noise could 
be a mode of vibration of a suspension system that couples in a suitable way the vibrations 
of mirrors in two different interferometers. It may be a valid assumption that all such sources 
of noise will be negligible. In this case we can use in Eq. ( |A33| ) the values = and 
Cmax = 1- [The value chosen for Cmax is unimportant as long as Cmax ^ 1/Vk]. Then from 
Eq. ( |A36| ) we see that (3abi ^ CmaxV^ as -C lO^®, and so to a good approximation. 
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pc{n) oc n n 

1=1 a<b 



- erf(-/5ab/) 



(A38) 



Essentially this pdf gives an upper bound on each Qj = Qg{fj) of fi/,max = miiia.b ^abi, 
which is of the order of ~ lO'^e if ^ 10^^. As already mentioned, this upper bound 
will be much worse than that obtained from intersite correlations, unless the dimensionless 
correlation coefficient e is ^ 10~^. Moreover, the bound is only obtained by making the 
specific assumptions (|A30|) and (|A37|) about sources of correlated noise at one site. For these 



reasons, in the following subsection on data-analysis we consider only intersite correlations 
and take pc{^) ~ const. 

We now show in more detail that very little information about the SB is obtained if the 
assumptions ( |A30| ) and ( |A371 ) are dropped. As a simple model, we consider the opposite 
extreme of assuming equal probability for all relative phases in the contributions to the 
outputs of different detectors from correlated sources of noise. Thus, we take the a priori 



pdf to be of the form |36 



PcoAc)dc = n n 

1=1 a<b 



X exp 



-abl\ 



2C2 



d(ReCabi)d(lmCabi), 



(A39) 



where Cmax is the a priori maximum correlation. Then from Eq. ( |A27| ) we obtain 



1=1 a<b 



(Re/?, 



abl 



1 + kC^^^ 

This function is of the form of Eq. ( [A42| ) below, where, from Eq. ( |A36|) , 



(A40) 



a, -10 W . 



(A41) 



This is roughly the minimum detectable value of Qg that can be detected in a bandwidth 
of ~ I/tc- Thus, the function ( |A40| ) is qualitatively similar to the pdf ( |A42| ) obtained 
from intersite correlations, with the simple change that the minimum detectable value of 



^Ig in any frequency band is increased by a factor of + kC^^^. If we take Cmax to be of 
order unity and so make no assumption about the amount of correlated noise present, then 
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the upper bounds on from intrasite correlations will be worse by a factor of \^ ~ 300 
than those obtained from intersite correlations. Only if Cmax ~ 1/ Vk will the two be 
comparable. However, because of the possibility of weak, unknown sources of noise, it is 
probably inappropriate to make such a strong assumption. 



4. Implications for data processing 



The distribution ( [A21| ) with pc = const is of the form 

p(^],)=A^3pf (^],) expJ 



1=1 



2a] 



(A42) 



where fli = ^g{fi = I/tc), because the matrices are all approximately diagonal in frequency. 
Using Eq. ( p.l2|) , we find that the argument of the exponential in Eq. ( |A21|) becomes 



^/^°"rf/tr{[S,-^ll-(n7^-S^)]'| 



(A43) 



where 7afe = 4pc7ab/ (5vr/^) . Hence we find, using df -* (I/tc) 'Zi, that 

1 



ktr 



(sl-'(//)-7)' 



(A44) 



and 



(A45) 



One would like to summarize the information contained in Eq. ( |A42|) by calculating 
some kind of signal to noise ratio. There are various, inequivalent ways of doing this. For 
example, one could calculate probability distributions for the quantities fimax = max/ Qj or 
Q = "^Zi^i/^'i it is clear that the SNR for fimax would be worse than that for Vt. What 
we shall in fact do is calculate the probability distribution p{VL\VLi = Vt) oi assuming 
that fii = . . . = f2„ = f2, as this is the easiest to calculate. Now if we ignore the a priori 
information represented by pf^ in Eq. ( |A42| ) [so that is a multivariate gaussian], then 

this is equivalent to calculating the pdf for the weighted average 
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(A46) 



However when we include the information contained in pf'\ the main effect is to truncate 
and renormahze |18[ the pdfs for each Qj, since all the f2's must be positive. Hence a simple 
interpretation of the pdf p{fl\Qj = Q) and the corresponding SNR in terms of the average 
( |A46| ) not possible. Nevertheless we suspect that p{Q\Qi = Q) approximately represents the 
probability distribution of some type of average of Qg{f). 

We now insert the assumption = Q (constant) into Eq. ( |A42|) , which gives a pdf of 
the form 



(A47) 



Here the quantity 



M 



I I I '^l 

is the statistic that should be calculated to estimate the value of Q. We obtain 



(A48) 



" -s 1 II 1 1 

OC / dftrSj^- " ■ 7^ ■ S^; 

Jo L 



(A49) 



This can be written in the time domain as 



^A/CX Jdt I dTHait + T)Labir)Hkit), 



(A50) 



where L(/) = 7(/)"'", and 



H(/) 



1 



fS/2 



(A51) 



The functions Labir) are the sliding delay functions discussed in Sec. |T|. Note that ^Im is 
constructed from the measured correlations in the following way: the intrasite correlations 
are used only to estimate the detector noise matrix S|, and then this noise matrix is used 
together with the measured inter-site correlations to estimate ^Im- However, the matrix S[[ 
will typically be approximately diagonal, and so neglecting the intrasite correlations will 
give only a negligible error. 
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Finally the intersite SNR VLm/c can be obtained from Eq. ( |A43|) ||18|| . If we assume that 
S/i(/)"'" = fircai(/)7(/)"'"5 so that all the intersite correlations we measure come from the SB, 
then the SNR takes the form 



1P=^{-^) Jo ^^—^ 



xtr 



(7(/)^-S;:^(/)")' 



(A52) 



Using the fact that the matrix '~f{f)ab will be constant on each subblock corresponding to 2 
detector sites, one can derive Eqs. ( p.l2| ) and ( |3.13| ) from Eq. ( [A52| ). 



APPENDIX B: CALCULATION OF THE 
OVERLAP REDUCTION FUNCTION 

The ath detector is characterized by its position and by the tensor d = (u (g) u — 
V ® v)/2, where u and v are unit vectors in the direction of its arms. In terms of these 
quantities, the overlap reduction function is given, from Eq. (p.l3|) , by 

7a6(/) = ^ E / (da : e^'") (d, : e^'") 

X exp [27rz/n • (x„ - x^)] . (Bl) 

If we write x^ — x;, = c? m where m is a unit vector, and put a = 2'Kfd (in units in which 
c = 1), then we obtain 

labU) = daij Tijki{a,m) dbku (B2) 

where 

r.,.K«, ^) = i-T.fd"^ e'--. (B3) 

This integral can be evaluated by the standard method of writing down the most general 
possible answer: 
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Tijkii a , m) = A{a)5ijSM + B{a) [SikSji + SuSjk] 

+ C{a) [dijrrikmi + dkimirrij] + D{a)mimjmk'mi 

+ E{a) [SikTrijmi + . . . + Sjiniimk] . (B4) 

One might expect to have to include a term proportional to wuWjk + WikWji where Wij = 
EijkTn^ 1 but in fact this tensor is a linear combination of the five tensors included above. 
Contracting Eq. (B^) with the five different tensorial expressions that appear on its right 



hand side yields a system of linear equations for A{a) . . . E{a) which involves scalar integrals 
that are straightforward to evaluate. Solving these equations and substituting the results 
back into Eqs. (|B2D and (|B4D gives 



labif) = Pi(tt)da : d;, + P2(tt)m ■ d„ ■ df, ■ m 

+P3(a)(m • da ■ m)(m ■ d;, ■ m). (B5) 

The functions pj {a) are linear combinations of the spherical bessel functions: 

pi(a) = 5jo(«) - 2ji(a)/a + 5j2(«)/a^ (B6) 



P2(a) = -lOjo(a) + 40ji(a)/a - 50j2(«)/a', (B7) 
and 

P3(«) = 5jo(«)/2 - 25ji(a)/a + 175j2(a)/(2a'). (B8) 

The result (^^) applies to any gravitational wave antennas, such as interferometers with 
non-perpendicular arms and arbitrary orientations, or resonant bar antennas. The special 
case of two resonant bar antennas, for which each d^ oc 1 — Sria for some vector n^, 
has been previously derived by Michelson []TB[. As can be seen from Fig. |TD|, the first term 
in Eq. ( [B5D dominates for a ^ 1 (unless d^ : d;, 0). A similar simplification applies for 
a ^ 1: we have 

latif) = 5jo(«) d,^ : d^ + 0{a-'), (B9) 
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where denotes the tracefree part of the projection {Sik — mimk){Sji — mjmi)daki of (!„ 
perpendicular to m. The fact that 7(/) does not depend on the components of and d^ 
parallel to m is due to the fact that the cross correlation at a ^ 1 is dominated by modes 
whose wavevectors are nearly parallel to m, as discussed in Sec. p.V B. 



To apply Eq. (BH) to terrestrial detectors, we need to choose a coordinate system and 



express the tensors di, d2 and m in terms of the angles 5, A and (3 defined in Sec. |V. 
A convenient choice is to take the detectors to be located at 9 = it/2 ± (3/2 and = 0, 
where r, 6, are spherical polar coordinates with origin at the earth's center, so that the 
unit vector in the direction joining the detectors is m = e^. If we let e^, and be the 
usual basis of orthonormal vectors, and define 



d(a, 9, 0) = sin(2CT)(eg ® - O e^)/2 

- cos(2a) (e^ ® + ® e^)/2, (BIO) 

then we can take di = d(A + 6,71/2 + f3/2, 0) and da = d(A -6,tt/2- (3/2, 0). Inserting 
these expressions into the formula ( |B5D yields, after some manipulation, the result 
where the functions gj{a) are 

^i(«) = ^/(«)-(-9,-6,9,3,l), (Bll) 

92{a) = ■ (45, 6, -45, 9, 3), (B12) 

lb 

g;{a) = ^f{a) ■ (15, -4, -15, 9,-1), (B13) 

and 

/(a) = (a cos a, cos a, sin a, o?' sin a, sin a) / oC' . (B14) 
From Eq. ( [4.1|) it is straightforward to evaluate the Fourier transform of the overlap 



reduction function, which gives the sliding delay function discussed in Sec. pi , 

/oo 
dfe-'-'^^jabif). (B15) 
-00 
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The function Lab{T) vanishes for \t\ > d = 2rEsm{/3/2) 
Eqs. (|4.1| ) to ([4.3D, where the functions gj are now 

and V = 1 — t'^/(P. 



For |r| < d, Labij) is given by 

(B16) 
(B17) 

(B18) 
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FIGURES 

FIG. 1. The broadband signal to noise ratio for a pair of detectors as a function of the angle (5 
subtended between them at the center of the earth, normalized to unity for coincident detectors. 
Curve I corresponds to each detector having an arm along the arc of the great circle that joins 
them, and curve II corresponds to each one having an arm at 45° to this arc. The optimal 
configuration is II for (3 ^ 70° (except very close to (5 = 0), and I for larger values of f3. The 
point C/L and the two points marked C/V show the expected sensitivities of the LIGO detector 
pair and of both LIGO/VIRGO detector pairs, with their current orientations. The point V/Q 
shows the sensitivity the detector pair VIRGO/GEO would have if the orientations were chosen 
optimally for the stochastic background (which will probably not be the case, since optimization 
for the stochastic background implies sensitivity to only one of the waves' two polarizations, and 
a corresponding loss of information when studying non-stochastic waves). 

FIG. 2. The "sliding delay function" L{t) for the LIGO pair of interferometers. To maximize 
the broadband sensitivity to the stochastic background, the cross correlation with a time delay r 
between the detector pair must be integrated against this function; see text. 

FIG. 3. The angles a\, (72, and P2 formed by a pair of interferometric detectors and the line 
L which joins them. 

FIG. 4. For coplanar detectors the overlap reduction function is of the form 
l{f) = cos(4(5)g'i(a) + cos(4A) [52(«) + fl'3(«)]; where 5 and A are the relative and total rota- 
tion angles, and a is the phase lag between the detectors at frequency /; see text. Here we plot 
the functions gi{a) and 5^2(0:) + gz{oL). 

FIG. 5. The overlap reduction function for the LIGO detectors using their currently planned 
orientations. 
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FIG. 6. For any orientation angles 5, A of the LIGO detectors, the overlap reduc- 
tion function will be given in terms of two functions of frequency Qi{f) and 62(/) by 
7(/) = cos(4(5)0i(/) + cos(4A)92(/). Here we show the functions 0i(/) and 02(/). 

FIG. 7. The broadband signal to noise ratio for the VIRGO/GEO detector pair as a function 
of the relative rotation angle 5, assuming (i) that the total rotation angle A is optimally chosen, 
and (ii) the advanced detector noise curve of Ref. The normalization is to unity for coincident, 
aligned detectors. 

FIG. 8. The functions 0i(/) and 02(/) appropriate for the LIGO detector in Hanford, Wash- 
ington together with the VIRGO detector in Pisa, Italy. The fact that 02(/) is large compared 
to 0i(/) for most frequencies implies that the sensitivity of the detector pair to the SB depends 
strongly in this case on the total rotation angle A of the two detectors with respect to the line 
joining them; see text. 

FIG. 9. A plot of the quantity f ^^Snif) "^ with arbitrary normalization as a function of fre- 
quency. It is this quantity that must be integrated against the square of the overlap reduction 
function to determine the broadband response of a detector pair to the stochastic background, 
when one assumes a constant, Harrison-Zel'dovich spectrum; see Eq. ( |1.2| ). It is clear that the 
highest sensitivity is limited to a narrow band of frequency between ~ 20 Hz and ~ 70 Hz. 

FIG. 10. A graph of the functions pi{a), P2{ct) and /93(a). 
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FIG. 1. The broadband signal to noise ratio for a pair of detectors as a 
function of the angle j3 subtended between them at the center of the earth, 
normalized to unity for coincident detectors. Curve I corresponds to each 
detector having an arm along the arc of the great circle that joins them, 
and curve II corresponds to each one having an arm at 45° to this arc. The 
optimal configuration is II for j3 < 70° (except very close to /3 = 0), and I for 
larger values of j3. The point Cj C and the two points marked CjV show the 
expected sensitivities of the LIGO detector pair and of both LIGO/VIRGO 
detector pairs, with their current orientations. The point V jQ shows the 
sensitivity the detector pair VIRGO^GEO would have if the orientations 
were chosen optimally for the stochastic background (which will probably 
not be the case, since optimization for the stochastic background implies 
sensitivity to only one of the waves' two polarizations, and a corresponding 
loss of information when studying non-stochastic waves). 




FIG. 2. The "sliding delay function" L[t) for the LIGO pair of interferome- 
ters. To maximize the broadband sensitivity to the stochastic background, 
the cross correlation with a time delay r between the detector pair must 
be integrated against this function; see text. 
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FIG. 3. The angles o"i, o"2, /3i, and j32 formed by a pair of interferometric 
detectors and the hne L which joins them. 
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FIG. 4. For coplanar detectors the overlap reduction function is of the form 
7(/) = cos(4^)5'i(q;)+cos(4A) [52(0=) + 53(0;)], where 6 and A are the relative 
and total rotation angles, and a is the phase lag between the detectors at 
frequency /; see text. Here we plot the functions 51(0;) and 52(0=) + 53(0;). 
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FIG. 5. The overlap reduction function for the LIGO detectors using their 
currently planned orientations. 
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FIG. 6. For any orientation angles 8, A of the LIGO detectors, the overlap 
reduction function will be given in terms of two functions of frequency 
0i(/) and 02(/) by 7(/) = cos(4^)0i(/) + cos(4A)02(/). Here we show 
the functions 0i(/) and 02(/). 
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FIG. 7. The broadband signal to noise ratio for the VIRGO/GEO detector 
pair as a function of the relative rotation angle 6, assuming (i) that the 
total rotation angle A is optimally chosen, and (ii) the advanced detector 
noise curve of Ref. [1]. The normalization is to unity for coincident, aligned 
detectors. 
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FIG. 8. The functions 0i(/) and 02(/) appropriate for the LIGO detector 
in Hanford, Washington together with the VIRGO detector in Pisa, Italy. 
The fact that 02(/) is large compared to 0i(/) for most frequencies implies 
that the sensitivity of the detector pair to the SB depends strongly in this 
case on the total rotation angle A of the two detectors with respect to the 
line joining them; see text. 
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FIG. 9. A plot of the quantity Sn{f)~^ with arbitrary normahzation as 
a function of frequency. It is this quantity that must be integrated against 
the square of the overlap reduction function to determine the broadband re- 
sponse of a detector pair to the stochastic background, when one assumes a 
constant, Harrison- Zeldovich spectrum. It is clear that the highest sensitiv- 
ity is limited to a narrow band of frequency between ~ 20 Hz and ~ 70 Hz. 




FIG. 10. A graph of the functions pi{a), P2{ct) and 103(0;). 
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